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Abstract 

We present high resolution measurements of the specific heat and the thermal 
expansion of the inorganic spin~Peierls cuprate CuGeOa in a magnetic field 
of 16 Tesla. At the transition from the incommensurate to the uniform phase 
both quantities show pronounced anomalies, which allow to derive the uni- 
axial pressure dependencies of the transition temperature. In high magnetic 
fields the specific heat is dominated by magnetic excitations and follows a 
T"^ law at low temperatures. The thermal expansion measurements show the 
occurrence of spontaneous strains along all three lattice constants and yield 
high resolution measurements of the temperature dependence of the incom- 
mensurate structural distortion. The sizes of the spontaneous strains in the 
incommensurate phase are significantly reduced, but both their anisotropy as 

well as their temperature dependencies are very similar to those in zero field. 
PACS: 65.70.+y,64.70.Kb,75.50Ee 
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During the last three years CuGeOs has been the subject of intensive investigations and 
the occurrence of a spin-Peierls transition (Tsp — 14K) in this inorganic compound is well 
established now. Most of the characteristic features of this transition are observed, e.g. the 
opening of a gap in the magnetic excitation spectrum below Tsp ^^-^ and the dimerization 
of the one-dimensional spin-i-Heisenberg chains In addition to the dimerization, 

very pronounced spontaneous strains (e) of all three lattice constants are found, which are 
proportional to the square of the structural order parameter Q [^0. This proportionality 
(e = kQ^) is expected from a Ginzburg-Landau expansion of the free energy when a usual 
linear quadratic coupling between Q and e is assumed (see e.g. refs. The temperature 

magnetic field phase diagram is also found to be in fair agreement with the theoretical 
expectations |[T0|JTT[1 . For low fields a reduction of Tsp occurs which is proportional to H^. 
For higher fields (H^12T) a transition to an incommensurate phase takes place. In this I 
phase the lattice distortion is expected to be incommensurate with respect to the underlying 
lattice. There are, however, different models for the spatial character of the incommensurate 
modulation; domain-walls or a sinusoidal modulation |]rT|. An incommensurate lattice 



modulation is indeed observed in CuGeOs for H^J 12T by x-ray scattering |TB|,|T^. Although 
these measurements favor the domain-wall picture, significant discrepancies between theory 



and experiment are still present ||14 



In order to study the magnetic and structural properties of the I phase we have performed 
high resolution measurements of such thermodynamic properties as specific heat and thermal 
expansion. Both quantities were measured on the same single crystal of CuGeOa with a size 
of about 6 X 5 X 8.3 mm^. The specific heat (Cp) was determined by a quasiadiabatic 
step-by-step method. The longitudinal thermal expansion coefficients = (l/Li)9Li/9T 
(Li denotes the length of the sample along the axis i) were measured with a capacitance 
dilatometer. 

Fig. ^displays the specific heat data of CuGeOs in a magnetic field of 16 Tesla which was 
oriented along the 6-axis. For comparison, we also show the data relative to zero magnetic 
field (see also refs. [p!5i|l6| ). It is obvious that the main infiuence of the field is a strong de- 



crease of both the transition temperature and the size of the anomaly. Nevertheless, a sharp 
and rather large anomaly of Cp is found at the transition between the I and the high tem- 
perature uniform (U) phases. The A-like shapes of the anomalies indicate that in CuGeOa 
the transitions are strongly affected by fluctuations, whereas the specific heat anomalies 
of the organic spin-Peierls compounds seem to show "mean-field" behavior [|T^-|TU[, i.e the 
"upturns" of Cp close to Tsp are not observed. There are two possibilities to discuss this dif- 
ference of the anomaly shapes. Either the better sample homogeneity of CuGeOs allows to 
observe the fluctuations at the spin-Peierls transition or fluctuations are more pronounced 
in CuGeOs. For example, one may speculate, that the good agreement with mean-field 
behavior in the organic compounds is related to the pre-existing soft phonon which has not 
been found in CuGeOs so far. 

Although the broadening of the transitions associated with sample inhomogeneities is 
rather small in CuGeOs, there is a significant asymmetry with regard to the maximum of 
Cp(T). Due to this asymmetry a description of the anomaly shape by a 'purely' critical 
behavior of Cp is questionable. As shown by the authors of ref. |]TB|] for their zero field data, 
which - although obtained on a different crystal - are in very good agreement with our 
data, it is not possible to unambiguously assign a universality class to the transition via the 
critical exponent of Cp. We note that - apart from their different sizes - the shapes of the 
anomalies close to Tsp do not differ significantly in H = and 16T. Thus, we do not analyze 
the critical behavior of Cp in this communication. 

In the low temperature range, however, we find distinct differences between the zero and 
high field data. The Inset of Fig. |l| shows Cp/T versus T^ for H=0 and 16 Tesla, respectively. 
As expected the data in zero field exhibit a clear curvature down to the lowest temperatures 
in this representation. This is due to the activated behavior of the magnetic contribution of 
Cp. The solid line represents a fit in terms of Cp = /9phT^ + 6 ■ exp(— AE/T). Restricting 
the fit to T < 6 K we obtain a value /5ph = 0.3 mJ mole~^ for the lattice contribution. 
The magnetic contribution is given by 5 = 3.6 J mole~^ and the energy gap amounts 
to AE = 23 K, in fair agreement with the results reported in ref. . 



The high field data show a completely different low temperature behavior. The entire 
specific heat follows a pure T^ law. There is no indication for a gap in the magnetic excitation 
spectrum of the I phase. Assuming that /3ph does not depend on H the specific heat is 
described by Cp = (/5ph + Pmag) ■ T^. For /^mag we determine a value of 1.4 mJ mole~^ K~^, 
which is significantly larger than /?ph. Thus, in the low temperature range the specific heat 
of the I phase is dominated by the magnetic excitations, which are strongly suppressed in 
zero field due to the opening of the gap. 

The T^ dependence of the magnetic specific heat is present in a rather large temperature 
range in the I phase, i.e. up to 5.5 K ZTsp/2. Moreover, we also find the T^ law with the 
same coefficient /?mag for smaller magnetic fields down to 12.5 Tesla. Apparently, neither 
the temperature dependence nor the absolute value of Cp strongly changes as a function of 
the magnetic field in the I phase. A very strong field dependence of both the temperature 
dependence and the absolute value of Cp at low temperatures is, however, present in the 
D-phase as can be extracted already from a comparison of the data in T and 6 T shown 
in ref . |T^ . 



Let us shortly discuss the implications of the T^ law we observe for the magnetic specific 
heat of the I phase. The - to our knowledge - only existing theoretical calculation of Cp, 



which uses the soliton lattice solution for the I phase based on a mean-field Hamiltonian pO 
can not explain our experimental finding. Fujita and Machida find in their calculations a 
BCS behavior of Cp for both zero field as well as in the I phase, i.e. they obtain a finite but 
reduced gap in the I phase as well. Even if we assume a very small value of this gap, it is 
impossible to fit our data with the theoretically expected BCS temperature dependence. 

From the experimentally observed specific heat one can in principle extract the disper- 
sion relations of the magnetic excitations. Assuming a Bose-statistics suggested by the T^ 
law and a dispersion relation of the form u; oc k"^ the specific heat at low temperatures is 
proportional to T^^"^ where D denotes the dimension of the system. For three-dimensional 
excitations our finding implies linear dispersion relations in the I phase of CuGeOs as in 
usual antiferromagnets. For D = 1 (or 2) rather strange dispersion relations are obtained. 



which makes the assumption of low-dimensional magnetic excitations in the I phase rather 
unlikely. Further theoretical and experimental work is necessary to decide, whether there 
are really three dimensional excitations with oc k in the I phase or whether the T'^ law of 
Cp we observe in a large temperature and field range is due to an accidental superposition 
of different contributions. 

Our findings for the thermal expansion in the incommensurate phase markedly differ 
from those of the specific heat. In Fig. ^ we show the thermal expansion coefficients along 
the a, b, and c axes. Concerning the thermal expansion one has to distinguish between the 
structural and magnetic anisotropy. The magnetic anisotropy arises from the different values 
of the gyromagnetic ratios along the different crystal axes. In our experimental setup 
the magnetic field is always parallel to the crystal axis whose thermal expansion coefficient is 
measured. Therefore, the transition temperatures slightly differ in the measurements carried 
out along the different crystal axes (see Fig. and Table |ID. 

The thermal expansion data in H = are also shown in Fig. ^ For each lattice constant 
the zero and high field data were obtained during the same run, i.e. for exactly the same 
orientation of the crystal. The anomalies occurring at the I/U transition are strongly re- 
duced in size (by about a factor 2) compared to those found in zero field. Despite this strong 
decrease, the anomalies are still rather large and indicate that the U/I transition strongly 
depends on pressure (see below). Furthermore, the anomalies reveal the occurrence of spon- 
taneous strains within the I phase, which are also seen in the temperature dependencies of 
the lattice constants (Fig. The spontaneous strains are obtained by integration of the 
anomalous contribution 6ai = ai — ai,extr. below Tsp. The Oi^extr. are smooth polynomials 
(given by the dashed lines in Fig. ^) representing the extrapolation of the cti in the U phase 
to at T = K 0. 

In zero field the €[ are proportional to the square of the order parameter P,[7|,p2| . The 



characteristic temperature dependence of an order parameter is also visible in H = 16 T 
(right part of Fig. The temperature dependence of the structural order parameter in the 
I phase, i.e. the incommensurate lattice deformation, has not been measured by diffraction 



techniques so far. Thus, we can not unambiguously prove ei oc in the I phase. However, 
this leading order of the strain order parameter coupling is usually observed at structural 
transitions and, moreover, there is no indication that higher order terms are important at 
the U/I transition in CuGeOs. In the following we will assume e oc for both phases and 
a comparison of the ei in H = 16T and H = allows for a comparison of the structural order 
parameters in the I and D phases, respectively. 

Besides differences in the critical behavior of the order parameter close to Tsp , which will 
be discussed elsewhere , the ei occurring at the I phase compare well to those in zero field. 
In Fig. ^ we show Tsp(H) ■ 6ai = — ^ versus reduced temperature t = 1 — T/Tsp(H). That 
means, we compare the derivatives ^ = k^^^^^ of the spontaneous strains occurring in the 
D and I phase, respectively, (k^''^ denote the strain order parameter coupling constants.) 
Please note that by comparing the temperature derivatives of the ei(t) possible differences of 
their temperature dependencies in H = and H = 16 Tesla will show up more clearly than 
by comparing the ei(t) themselves. It is apparent from Fig. ^that, apart from the absolute 
values (see the different scales in Fig. the temperature derivatives of the in the I phase 
are very similar to those in H = 0. This holds for all three crystal axes. Moreover, the ratio 
between ^ in H = and in H = 16 T is independent on the crystal axis and amounts to a 
value of 3, i.e. the structural anisotropy does not change in a magnetic field. We emphasize 
that for the specific heat in H = and H = 16T such a similarity is not present at all 
(see Fig. [l|). The magnetic and the structural degrees of freedom show completely different 
magnetic field dependencies. 

The reduced ei in the I phase may arise either from three wave vector dependent strain 
order parameter coupling constants (k[ = kj(q)) or a reduced average amplitude A of the 
distortion. The identical anisotropy in H = and H = 16T strongly indicates that the 
reduction for all three lattice constants is determined by a common parameter, e.g. by A. 
There is a straightforward qualitative explanation of such a reduction within the domain 
wall picture. A is expected to reduce when approaching the domain walls [O . The reduced 



ei in the I phase is then a consequence of the reduced average value of the dimerization 
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within a single domain. Moreover, the ei should show a further decrease with increasing 
field due to the increasing number of domains. However, a sinusoidal modulation of the 
lattice distortion (A(x) = Aosin(qx)) also causes a reduction of the average value of the 
dimerization. Thus, from our present data it is impossible to discriminate between these 
two models of the structural distortion. Measurements up to still higher fields are planned 
to study the field dependence of ei(T, H) in the I phase to clarify this question. 

Finally we derive the uniaxial pressure dependencies of the U/I transition by comparing 
the anomaly sizes of the ai and Cp ||23| as described in ref. [|7|,p2|. Similar as for H = T the 



uniaxial pressure dependencies of Tsp in H = 16T are very large and strongly anisotropic (see 
Table D . For pressure along the a axis Tsp decreases, whereas it increases for pressure along 
the two other axes. The hydrostatic pressure dependence given by the sum of the 9Tsp/9pi 
amounts to +3.6K/GPa. All the pressure dependencies in H = 16T are significantly smaller 
than those in H = (Table j!]). There are, however, pronounced similarities to the zero 
field data. The anisotropics of the dTsp/dpi at the U/I transition and at H = are nearly 
identical. Moreover, the pressure induced relative changes of Tgp in H = and H = 16T 
are very similar. For instance, along the b axis an increase of about 50% /GPa is found for 
both Tsp(H = 0) and Tsp(H = 16T) indicating a pressure independent H/T phase diagram 
in reduced scales, which is in agreement with theoretical predictions. However, as we have 



recently shown p4| the (9Tsp/i9pi unexpectedly correlate with the magnetoelastic coupling. 

To summarize, we have reported measurements of the specific heat and the thermal ex- 
pansion of CuGeOs in H = 16T. For both quantities we find pronounced anomalies at the 
U/I transition. The specific heat at low temperatures is dominated by the magnetic excita- 
tions. In contrast to the findings in zero field our data in H = 16T yield no indication for a 
gap in the magnetic excitations. Instead the magnetic specific heat follows a pure T^ law. 
The thermal expansion below Tsp is dominated by the incommensurate lattice distortion 
leading to spontaneous strains of all three lattice constants which scale with the structural 
order parameter. Remarkably, the temperature dependencies of the commensurate and in- 
commensurate lattice distortion at H = and 16 T, respectively, are very similar, whereas 



the magnetic specific heats markedly differ. 

The work at Cologne was supported by the Deutsche Forschungsgemeinschaft through 
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FIGURES 

FIG. 1. Specific heat of CuGeOs in H=0 (o) and 16 Tesla (•). Inset: Cp/T vs. T^; the 
sohd hnes are fits assuming activated behavior (o, H = 0) and a T"^ law (•, H = 16T, see text), 
respectively. 

FIG. 2. Thermal expansion coefficients of CuGeOa in H=0 (•) and 16 Tesla (o). The dashed 
lines represent the extrapolated behavior of cci of the U phase to ai{T = 0) = 0. The arrows mark 
the transition temperatures derived at the maximum slopes of ai{T) 

FIG. 3. Left panel: Temperature dependence of the lattice constants of CuGeOs in H=16 Tesla. 
Right panel: Spontaneous strains ei = / ai — ^ai^extr. in the I phase. 

FIG. 4. Temperature derivatives of the spontaneous strains — ^ oc vs. reduced tempera- 
ture t in H = 16T (o, left scale) and H = (•, right scale). 
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TABLES 



axis 


H=0 Tsp=14.35(10)K 


H = 16 Tesla 




^ |p,=o (K/GPa) 


Tsp (K) 


^ |p,=o (K/GPa) 


a 


-3.7(5) 


10.13(10) 


-2.8(5) 


b 


7.2(5) 


10.02(10) 


5.3(5) 


c 


1.6(5) 


10.23(10) 


1.1(5) 



TABLE L Uniaxial pressure dependencies of Tgp along the different crystal axes in H = and 
H = 16T as calculated from the thermal expansion and specific heat anomalies The different 
Tgp's in H = 16T are due to the slightly different g-factors of CuGeOs. 
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